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Abstract 

Momentum production during the phase shp process in SNS Joseph- 
son junction is discussed. It is caused by the spectral flow of bound 
states of fermions localized within the junction. This effectively re- 
duces the Magnus force acting on vortices which provides an explana- 
tion for the experimental observation of the negligible Magnus force 
in 2D Josephson junction arrays. The flow of the fermionic levels is 
similar to that in sphalerons in particle physics, where it gives rise to 
the baryogenesis. 
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Vortex dynamics in 2D Josephson junction arrays has been studied in- 
tensively in recent years. The vortices in the arrays could be considered as 
massive particles with long-range Coulomb interaction 0]. A particular 
attention was devoted to the forces acting on the vortices, one of them being 
Magnus force. 

In the experiment the straightforward ballistic motion of vortices was 
observed which implies the absence of reactive forces acting on a vortex 
perpendicular to its velocity. This was also confirmed by more recent exper- 
iments: vortices move perpendicular to the driving current @], and no Hall 
effect was detected in the system p. 

An explanation of the absence of the Magnus force was proposed in 
Ref. 0. The authors claim that the Magnus force is proportional to the 
offset charges on the superconducting islands, the effect of which being neg- 
ligible. In a recent worked [[7[ Gaitan and Shenoy argued that the Magnus 
force is proportional to the density of superconducting electrons on the is- 
lands averaged over distances large compared to the lattice constant of the 
array rather than to the charge of the island which is given by the difference 
in the numbers of electrons and protons. On the other hand Zhu, Tan and 
Ao have shown the force to be proportional to the local superconduct- 
ing density at the point where the vortex is situated. Since the vortex does 
not move through superconducting islands but through the junctions, the 
Magnus force on the vortex can be substantially reduced. 

In the present paper we consider the forces on a vortex moving in a 
Josephson junction array, and propose a different explanation to the exper- 
iments mentioned above. It should be emphasized that one contribution to 
the force was missed in previous considerations 0, ^ §] which is the force 
from the spectral flow. In uniform superfluids there exists the so called zero 
branch of energy levels of fermions localized within vortex core which crosses 
zero as a function of angular momentum . Vortex motion leads to the flow 
of the levels along this branch, and energy of some levels crosses zero value. 
At low temperature these levels become occupied (or unoccupied, depending 
on whether they cross zero upwards or downwards). During this process the 
whole number of quasiparticles localized within the core is conserved, while 
the linear momentum of the quasiparticles is not conserved. This implies a 
transfer of the linear momentum from the vortex to the heat bath and thus 



an additional force acting on the moving vortex |10[. This force from the 



spectral flow can almost completely neutralize Magnus force [O] though in 
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some regimes the spectral flow is suppressed and the net force appears. This 
scenario reproduces the microscopic calculations by Kopnin and co-authors 

[0- 

We argue that an analogous situation takes place in Josephson junction 
arrays. Here the role of fermions localized within a vortex core is played 
by fermions localized within junctions. The phase slip events in junctions 
during vortex motion are accompanied by the spectral flow of the bound 
state fermions. The force from this spectral flow neutralizes Magnus force, 
and the reactive force on a vortex becomes negligible. 

We consider a particular example of a Josephson junction array and show 
the cancelation of two contributions to reactive force. Other geometries and 
possible difference in Magnus and spectral flow forces are under investigation. 

We consider a square lattice of superconducting islands as plaquettes, the 
edges of the lattice representing normal layers in SNS-junctions. Vortices 
could be considered as sitting at vertices of the lattice (Fig.]!]). 

As in the review [|l3l and Ref.[jlll (see also fl^) we assume that the 



conventional Magnus force on a vortex at low temperature is determined 
by the local density n of the electron because n is the variable which is 
canonically conjugated to the phase of the Bose-condensate: 



' M 



Tcnz X , (1) 



Here is the velocity of the vortex with respect to the superfluid velocity, 
which is chosen to coincide with the velocity of the heat bath (v^ = v„ = 0). 
The 2D particle density of electrons n is assumed to be the same in the 
superconducting and normal regions. Note that in the low-temperature limit 
the superfluid density Ug tends to n irrespective of the magnitude of the order 
parameter gap A (actually due to the lordanskii force from the heat bath, 
the Eq. ([|) is valid even at nonzero T [0, |1^ ) . This can be applied also to the 



electrons in "normal" regions with small magnitude of A. This expression 
does not contradict to Ref.[0] at T = 0, but we want to demonstrate that 
there is another force which nearly completely compensates the Magnus force. 

During one step a vortex moves from a vertex to an adjacent one. We are 
going to calculate production of momentum due to spectral flow of fermions 
for such a process, which gives rise to the compensating force. The problem is 
in many details similar to the evolution of fermionic bound states, which exist 
within topological solitons in polymers, superfluid ^He and other ordered 
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systems in condensed matter with fermions |T5[. The same type of evolution 



of the fermionic spectrum in topological solitons and sphalerons in particle 
physics leads to baryogenesis JTB], , while in our case the spectral flow leads 
to " momentogenesis" . 

Let us consider one particular Josephson junction with x being the co- 
ordinate axis normal to the junction (Fig.|l]). The dependence of the eigen- 
functions of Bogolyubov-Nambu hamiltonian 

'H = %-e + fiReA - fglmA (2) 

on r_|_ = (y, z) is given by a factor exp(2k_|_r_|_). Here r are Pauli matrices in 
Bogolyubov-Nambu space and e = (— — kp)/2m* is the operator of energy 
of quasiparticles in normal liquid; kp is the Fermi momentum. We suppose 
that the order parameter varies slow on the length scale of the coherence 
length, while the Fermi momentum remains constant. We are interested in 
the energies close to Fermi surface, so, we suppose that k±_ < kp and the 
eigenfunction is represented in eikonal approximation: 

exp{iqx) ■ ip{x) (3) 

where = kp — k'j_ {q plays a part of Fermi momentum of ID Fermi liquid for 
given k±), the exponent represents fast oscillations in space, and tl^lx) varies 
slowly. In this case we may substitute e by q{—iV)/m*. The approximation 
gives energy spectrum with the accuracy small compared to A. 

We shall investigate the dependence of the energy spectrum on the phase 
difference A0 between two superconducting islands. For simplicity we sup- 
pose that the order parameter is real A(x) = |A(x)| on the one side of the 
junction {x < 0), and on the other side [x > 0) the order parameter is given 
by A(a;) = |A(x)| ■ exp(i0). For (p = tt the order parameter is an odd real 
function of x (Fig.^ and the hamiltonian is supersymmetric: {H, T2} = 0. 
This — TT soliton corresponds to the sphaleron in particle physics [T^, |l^], 
i.e. the intermediate state between vacua with different topological charges. 
Due to supersymmetry the Hamiltonian has an eigenfunction 

i>Q = const ( . / V Wo, 
\^ -sign(q)z J 

V'o = exp I -^|da;'A(x') I . (4) 
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with zero eigenvalue. So, at = vr for each q one energy level crosses zero. 
For small — tt the perturbation of the hamiltonian is 



H^n^ = -(0-vr)A(x)f2 (5) 
for a; > and Tiint = for a; < 0. The energy level is shifted to 

^ = sign(g)(0-7r)a;(|g|), (6) 

where 

oo 

S dx\A{x)\xljl{x) 
^(kl) = '-^ • (7) 



From it follows that at = vr, i.e. at the sphaleron, the energy level 
crosses zero upwards for g > and downwards for g < 0. The similar 



phenomenon (cf. Refs.|]l8|, ^ as well) was found for sphalerons in particle 
physics. 

This leads to the production of the x-component of the linear momentum 

AP = 2--A/^2|g| = SA (8) 

The prefactor 2 stands for double spin degeneracy, and 1/2 removes the 
double counting of particle and hole momenta. Here A is the area of the 
cross-section of the junction, which is given by the product of the lattice 
constant a and the thickness h of the film along z-axis (which is perpendicular 
to the plane): A = ah. 

Using this equation one can find the spectral-flow force experienced by 
the vortex moving with respect to the heat bath. Since the velocity of the 
vortex is vl = a/ At, where At is the period of time during which the vortex 
crosses one junction, one has for the spectral-flow force in x direction 

AP APvL „ 

^sp.flow = = = TTGoI'l, (yj 

where Cq = kkp/Sn"^. In vector notations 

Fsp.flow = ttCqVl X z. (10) 
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The 2D particle density is n = ph where the 3D density p ~ kp/Sn^ with 
the accuracy of the order of {/S./Ep)'^ . Therefore, the force induced by the 
spectral flow nearly neutralizes the Magnus force as anticipated. 

Discussion. Let us now make some remarks on a more general case of 
the geometry of Josephson junction array when superconducting islands do 
not cover almost all area of a sample. If the core size of a vortex is small 
compared to the lattice spacing a then one has the situation discussed in 
the vortex can be considered as a point-like object moving in the locally 
homogeneous environment with the slowly changing density of the electrons. 



In this case one can apply the bulk results [|r^: the en route Magnus force 
will be canceled by the en route spectral-flow force for any route. 

Let us consider another extreme limit Tc 3> a, which takes place if su- 
perconducting islands cover a small part of the area. In this case considered 
in Ref. 1^ one can average over the distances of the order of the core size 
and obtains the case of conventional homogeneous superfiuid with as a 
coherence length and a as "interatomic distance" . The Magnus force in this 
case will be determined by the average number density n as it was stated in 
Ref. 0, however the spectral flow force on a vortex will be also proportional 
to average value of Cq. Thus they again will nearly compensate each other, 
but now with the relative accuracy of the order of (a/rc)^. 

In the intermediate case of rc ~ a the difference in two contributions 
to the reactive force can become comparable to their values, because the 
en route values of n and Co can be different. Note that there are at least 3 
different reasons which can violate the fine tuning between n and Cq: (1) The 
particle hole asymmetry gives the disbalance n — Cq ^ n{A/ EpY . (2) The 
inhomogeneity discussed above gives n — Co ~ n{a/rcY. (3) Finite spacing 
ojq between the quasiparticle energy levels, which supresses spectral flow of 
fermions: n — Cq ~ niujoTY 0- In conventional Abrikosov vortices the 



spacing is c^o ~ ^^/Ep, while in the Josephson junction the discreteness 
arises due to the dimensional quantization along the axis y, which is parallel 
to the junction. 
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Figure 1: Josephson junction array. The squares represent superconducting 
islands and white spacings between them are normal intermediate layers in 
junctions. The distribution of phases of superconducting islands for a vortex 
sitting on a site of the lattice marked by a black circle is shown by arrows. 
When the vortex moves to the site marked by a dashed circle some levels of 
fermions localized within the junction between two central squares cross zero, 
and some amount of linear momentum kx is produced. This corresponds to 
an additional force on the vortex which compensates conventional Magnus 
force. 



Figure 2: The dependence of superconducting order parameter on the coordi- 
nate X normal to a junction for phase difference of tt between superconduct- 
ing islands. Corresponding Bogolyubov-Nambu hamiltonian has zero energy 
eigenvalues, i.e. at the moment when crosses tt the levels cross zero. 



